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PLANAR PARTICLE/DROPLET SIZE 
MEASUREMENT TECHNIQUE USING DIGITAL 
PARTICLE IMAGE VELOCIMETRY IMAGE DATA 

This invention was made with government support under Grant No. 
NGT 3-52352 NASA Glenn Research Center and Grant No. NAG3-2110 NASA 
Glenn Research Center, 

Background of the Invention 

5 The present invention relates to particle sizing. It finds particular 

application in conjunction with two-phase flows and will be described wdth particular 
reference thereto. It will be appreciated, however, that the invention is also amenable 
to other like applications. 

Two-phase flows are used in many fuel combustion processes such as 

10 those found in gas turbine combustors, coal furnaces for power generation, and diesel 
engines. The efficiency of two-phase flows directly impacts the efficiency of the 
associated combustion process. For many industrial and fossil fuel energy processes, 
both the droplet/particle size and spatial distribution are of interest. The useful 
diagnostic tools should be able to make in situ measurements without disturbing the 

1 5 flow field so that the measurements can be meaningful for understanding these flows. 
This suggests the need for a novel non-intrusive optical technique that can provide 
instantaneous measurements of particle size and velocity at multiple spatial points in 
planar (2-D) fields so that estimates of the mass flow are obtained. 

Digital Particle Image Velocimetry (DPIV) is a technique for obtaining 

20 planar measurements of particulate seeded flow fields. Light from a pulsed laser is 



formed into a thin sheet to illuminate a planar cross section of a flow. A CCD camera 
is used to record the light scattered by the particles at the two (2) instants that the light 
sheet is pulsed. The fluid velocity is determined by analyzing the recorded particle 
image data. Typically, cross-correlation data analysis is used to reduce the recorded 
5 image data to determine the fluid velocities. In most instances, the flow field is 
artificially seeded vnth tracer particles in order to measure the flow velocity. In two- 
phase flows, the second phase material provides the scattering sites for estimation of 
the flow velocity, or at least the velocity of the second phase. The objective of this 
work is to determine the feasibility of estimating particle size fi-om Particle Image 

10 Velocimetry (PIV) image data. Particle size information can most likely be extracted 
provided the imaged particle size exceeds the optical system blur circle. The other 
major factors affecting the accuracy to which particle size information can be extracted 
are the optical system fi^number, image system pixel resolution and dynamic range, the 
optical properties of the particles, and the characteristics of the scattered light fi"om the 

1 5 particles. 

Many different optical techniques for making in situ measurements of 
particle/droplet sizes are knovm. 

An in-line holography system has been used to record holographic 
images of the light scattered fi-om particles in a fluid. The hologram has been shown 
20 to be that of a screen with a circular aperture representing the particle image. 
Holograms have been reconstructed by passing a laser beam through the hologram and 
recording the reconstructed image at an on-axis observation plane. The reconstructed 
image contains the Fraunhofer pattem of the imaged particles. The size and shape of 
the particle is then determined by direct observation of the reconstructed image. 
25 Holography has been easily extended to velocity measurement by using double 
exposure holograms and analyzing the double aperture type of fringe structure in the 
reconstructed image to calculate the particle separation. 

Laser interferometry has been used to characterize particles using the 
spatial frequency of the far-field fringe pattem in the forward-scatter region. Ovryn 4 




uses partially coherent light and forward scattering to obtain Poisson Spot images of 
particles suspended in a solution to determine their 3-D velocities. Although not a 
sizing technique, this work illustrates the information content in the diffraction rings 
surrounding coherently illuminated particles and such rings are used to determine 
5 additional properties about the particle. The difficulty with these holographic and 
interferometric techniques lies in the requirement of a high spatial resolution detection 
media, typically photographic film. 

Phase Doppler Particle Analysis (PDPA) is an existing technique that is 
able to make point wise velocity and size measurements simultaneously. Droplet sizes 

10 are obtained using PDPA, which relates the droplet size to the phase shift of light 
refi-acted through the drop and scattered to different positions on the receiving lens. 
The technique known as Laser Doppler Anemometry (LDA), which is a technique for 
obtaining point wise velocity measurements, has been modified for applications to 
particle sizing as well. The particular modifications relate the particle diameter to the 

15 LDA signal visibility. Other variations of the technique relate either the scattered 
intensity or the phase shift of the scattered light to the particle diameter. 

In diffraction-based sizing techniques the particle interaction with the 
illumination is assumed to be analogous to the interaction of the same illumination 
v^th a imiformly illuminated circular aperture. The inherent problem in diffraction 

20 methods lies in the assumption that a particle acts as a uniformly illuminated aperture, 
which ignores many other scattering effects present in light scattering firom a spherical 
particle, such as specular reflection, interference and refiraction effects. Therefore, the 
Mie scattered signal has been used for determining particle characteristics. 

The technique of Laser Sheet Dropsizing (LSD) has been used to 

25 determine droplet diameters. LSD uses the Mie signal, which is proportional to the 
particle diameter squared, and the Laser Induced Fluorescence (LIF) signal, which is 
proportional to the particle diameter cubed, for determining particle size. 

Mie scattering has been investigated from small particles in the regime 
where the particle image is dominated by diffraction to determine the minimum 




particle diameter that is detectable for use in pulsed laser velocimetry techniques. 
Taking the Fourier transform of the Mie scattered electric fields from spherical 
particles results in multiple peaks in the intensity profile, the characteristics of which 
can be related to the particle size. 
5 The ratio of the projected area of the reflection glare spot to the cross- 

sectional area of the droplet has been found to be proportional to the square of the ratio 
of the aperture radius to the distance from the droplet to the sensor surface. The signal 
to noise ratio of the light scattered from a spherical water droplet of known size is used 
as a measure of the ratio of light scattered from the glare spot to light scattered by the 

10 other part of the sphere surface. This information is then used to predict the shape and 
size of the reflection glare spot. Hence, information contained only in the reflection 
glare spot has been related to the droplet size. However, until now particle sizes 
(diameters) have not been determined as a function of the distance between a 
reflection glare spot and a transmission glare spot. 

15 The present invention provides a new and improved apparatus and 

method which overcomes the above-referenced problems and others. 



Summary of the Invention 

A method for determining a mass flux of an entrained phase in a planar 
two-phase flow records images of particles in the two-phase flow. Respective sizes of 
20 the particles are determined as a function of a separation between spots identified on 
the particle images. Respective velocities of the particles are determined. The mass 
flux of the entrained phase is determined as a fimction of the size and velocity of the 
particles. 

In accordance with one aspect of the invention, the recording step 
25 includes recording an image of a transparent particle. 

In accordance with another aspect of the invention, glare spots are 
identified on the particle. The particle size is determined as a function of a separation 
between the glare spots. 
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In accordance with another aspect of the invention, the step of 
determining the velocity includes determining the velocity as a function of a 
velocimetry of the particles within the images. 

In accordance v^dth a more limited aspect of the invention, the step of 
5 determining the velocity as a function of the velocimetry includes obtaining two 
exposures of the respective glare spots and measuring a displacement between the two 
exposures during a specified time interval. Either the average displacement 
(correlation processing) or individual particle displacements can be estimated 
depending on the data processing technique applied and the concentration of seed 
10 psirticles. 

In accordance with another aspect of the invention, the step of 
determining the velocity as a function of the velocimetry includes detecting a Doppler 
shift of light. 

One advantage of the present invention is that it makes in situ 
15 measurements without disturbing a flow field. 

Another advantage of the present invention is that it provides a non- 
intrusive optical technique that offers instantaneous measurements of particle size and 
velocity at multiple spatial points in planar (2-D) fields, thereby permitting an estimate 
of mass flow to be obtained. 
20 Still further advantages of the present invention will become apparent 

to those of ordinary skill in the art upon reading and tmderstanding the following 
detailed description of the preferred embodiments. 

Brief Description of the Drawings 

The invention may take form in various components and arrangements 
25 of components, and in various steps and arrangements of steps. The drawings are only 
for purposes of illustrating a preferred embodiment and are not to be construed as 
limiting the invention. 

-5 - 




FIGURE 1 illustrates scattering angles, 6 and 0, used in a Mie 
Scattering algorithm; 

FIGURE 2 illustrates a scattered intensity as a function of 6 for a 100 
p.m diameter PSL sphere in water with 532 nm wavelength illumination; 
5 FIGURE 3 illustrates the production of transmission and reflection 

glare spots; 

FIGURE 4 illustrates glare spot separation as a function of particle 

diameter; 

FIGURE 5 illustrates an algorithm for scanning PIV images and 
10 determining the size of detected particles; and 

FIGURE 6 illustrates glare spot intensity fields for a range of particle 

diameters. 

Detailed Description of the Preferred Embodiments 

A technique making use of the characteristics of the glare spots 
15 predicted by Mie theory is used in the present invention for particle sizing using 
Digital Particle Image Velocimetry (DPIV). A DPIV system is used to image the light 
scattered from spherical particles under controlled conditions at an observation angle 
of 90^. This recorded signal is representative of the Fourier transform of the Mie 
scattered electric fields incident on the lens of the imaging system. An evaluation of 
20 the present invention has been performed via a computer simulation program, which 
incorporates Mie theory and the Fourier transform properties of lenses on the scattered 
electric field. The location of the glare spots resulting from reflected and transmitted 
rays in the direction of observation are an indicator of particle size. While this 
measurement technique is not as accurate as point measurement techniques such as 
25 Phase Doppler Particle Analysis (PDPA), it has the potential for providing reasonable 
accuracy for size and high accuracy for velocity estimates in planar 2-D fields so that 
mass flow rates can be determined. The determination of mass flow rates are essential 
in many production type facilities (e.g., coal furnaces). Experiments have been 



performed to verify the simulation and to determine the feasibiHty of the present 
technique for measuring particle size from Digital Particle Image Velocimetry (DPIV) 
imagery. 

Mie Theory Background 
5 Mie theory provides the exact solution of electromagnetic waves 

scattered from a spherical particle that includes all of the diffraction, specular 
reflection, interference, and refraction effects. At an observation angle of 90^, two (2) 
glare spots are distinguishable on the surface of a transparent absorbing particle. The 
positions of these glare spots are described in Mie theory as a function of the particle 
10 diameter, the relative refractive index, and the angle from the optical axis at which the 
particle is being viewed. Knowledge of these positions is used to characterize the 
particles. 

The Mie scattering problem defined here assumes plane wave 
— > 

illumination v^th the electric field E polarized in the x-direction (perpendicular to the 
15 scattering (y-z) plane) incident on a homogeneous, absorbing sphere in a 
non-absorbing medium. This polarization is referred to as s-polarization or transverse 
electric (TE) polarization. With reference to FIGURE 1, a particle 10 within a two- 
phase flow is assumed to be in focus at the CCD detector plane. Angles G and (j) are 
defined from the x and z axes to a scattering plane 12 in which a detector (e.g., CCD 
20 camera) 14 is located. Ein^, Bi„c> and lq„j. are the electric and magnetic fields and the 
propagation direction of the incident illimiination. The detector 14 communicates with 
a processing device 16, which, as discussed below, calculates the size and velocity of 
the particle 10. 

The Mie scattering solution has been derived for the problem of the 
25 scattering of plane electromagnetic waves from a sphere. The magnitude of the 
scattered intensity is defined in Equation 1 as: 
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where Eq is the magnitude of the incident electric field, S,, and are the scattering 
amplitudes for the ^ and 0 polarizations of the scattered light, respectively, is the 
permeability of fi"ee space, c is the speed of light, r is the distance fi-om the particle to 
5 the detector plane, and k is the propagation constant. This relation shows the 
polarization dependence of the scattered electric field firom a particle at the observation 
angles 0 and <j). In the case of PIV imaging, the G range is about 90** and the 9 and (f) 
range is about 0** or 90®, depending on whether the incident electric field is parallel 
(p-polarized, or transverse magnetic (TM) polarized incident light) or perpendicular 

10 (s-polarized or TE polarized incident light) to the scattering plane. When ^ = 0°, the S, 
contribution is negligible since sin 0° is equal to zero (0). Conversely, when <|) = 90° 
the S2 contribution is negligible. For s-polarized light = 90°), the reflection and 
transmission glare spots are approximately equal in amplitude and this is the preferred 
illumination state. For p-polarized light (([> = 0°), the transmission glare spot is 

15 significantly dimmer than the reflection glare spot, which is not the preferred 
condition for obtaining accurate estimates of the glare spot peak locations. FIGURE 2 
shows how the Sj and 83 electric fields 20, 22, respectively, vary over a typical range 
of 6 values. 

When the incident light is unpolarized, or randomly polarized, the 
20 magnitude of the scattered intensity is defined in Equation 2 as: 

which has equal contributions of S,, and 83 scattered electric fields. 

The light scattered by the particles is imaged onto a detector in order to 
make a measurement. The work done on the electric field imaged onto a CCD array 
25 by a lens is analogous to taking the Fourier transform of the input electric field. 
Taking the Fourier transform of the electric field signal and multiplying it by its 
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complex conjugate yields the light intensity distribution on the CCD array (see 
Equation 3). 

^scattered ~ ^ (^scattered )^ ^ (^scattered ) (^) 

where 3 is the Fourier Transform, * represents the complex conjugate, hcattered is the 
5 scattered intensity imaged on the CCD detector and Escmtered the scattered electric field 
at the lens plane of the optical system. This form of the signal makes the scattering 
effects, such as reflection and transmission, visible in the form of glare spot structures 
in the particle image. 

Glare Spots 

10 With reference to FIGURES, glare spots 30, 32 are bright spots 

resulting from off-axis light scattering by a spherical particle 34. According to ray 
theory a glare spot is seen at the exit points of any rays reflected or refiracted in the 
direction of the observer. FIGURE 3 demonstrates ray theory for the reflected rays 36 
and transmitted rays 38 at a 90^ observation angle. The points 42 and 44 correspond to 

15 the reflected and transmitted glare spots 30, 32, respectively. The ray denoted as the 
reflection ray 36 hits the surface of the spherical particle 34 and is reflected fi:om the 
surface in the direction of observation. The ray denoted as the transmission ray 38 is 
intemally reflected once and is then transmitted out of the particle 34 in the direction 
of observation. 

20 The distance Xq + Xj between the glare spots 30, 32 is directly related 

to the diameter of the particle 34. Therefore, it is possible to use the glare spot 
characteristics to determine particle/droplet sizes. As the particle diameter increases 
so does the separation distance between the glare spots. The glare spots 30, 32 are 
oriented in the direction of propagation of the incident light. The location of glare 

25 spots is calculated based on the size parameter (27ia/X, where a is the particle radius 
and X is the illumination wavelength), and the observation angle, Oq, by studying the 
Fourier transform of the geometric scattering amplitudes and the Mie scattered electric 



-9- 



fields. The essence of the invention consists of estimating the particle size from a 
particle image on a CCD detector. 

The glare spot positions relative to the particle centroid, X/^j, in the 
object plane are determined and then related to the size parameter and observation 
5 angle. Therefore, if the positions are known at a specific observation angle, the size of 
the particle may be determined. The positions of the reflection and transmission glare 
spots on the detector plane are defined as: 

X^ = -aM cos^ ... (4) 

(5) 

10 where n is the ratio of the mdex of refiraction of the particle material to the index of 
refi-action of the mediimi, a is the particle radius, and M is the optical system 
magnification. The measured separation between the glare spots on the CCD detector 
leads to determination of the particle diameter dp via the following formula: 
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15 where A is the number of pixels separating the glare spots, Sp is the CCD pixel size, M 
is the optical system magnification, n is the relative index of refi*action, and 0q is the 
observation angle. This relation is used to analyze data obtained from PIV images to 
estimate particle sizes in a flow field. 
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Using Glare Spot Separation To Estimate Particle Diameter 
20 A one-to-one correspondence exists between the particle size and the 

glare spot separation according to Mie theory. Equation (6) gives the particle diameter 
in micrometers as a function of glare spot separation, A, in pixels, CCD pixel size, 8p, 
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in micrometers per pixel, relative index of refraction, n, optical system magnification, 
M, and observation angle, G^. FIGURE 4 shows a plot 50 of glare spot separation as a 
function of particle diameter based on this equation for a polystyrene latex sphere in 
water ranging in diameter from 45 yim to 160 \xm with M = 0.31, Sp = 6.7 jom/pixel, 
5 n= 1.195, and 0o==9O**. The separation with respect to particle diameter is linear 
showing a one-to-one relationship between the glare spot separation and the particle 
diameter. Therefore, the glare spot separation is a viable indicator of the particle size. 
However, due to the low spatial resolution of a CCD camera, it is still difficult to 
extract the particle sizes to high accuracy based on the spot separation. 

10 The accuracy of the particle size estimates is determined by the 

accuracy by which the two (2) glare spot peaks are determined. It has been shown that 
the uncertainty in estimating the center location of a Gaussian signal is proportional to 
the ratio of the width of the signal to the square root of its amplitude. In standard PIV 
data processing, the particle images or correlation peaks are typically estimated within 

15 an accuracy of about 0.2 pixels, provided the particle images span 1-2 pixels across 
their width. The glare spots fi-om the particle images satisfied this criterion. Hence, 

the glare spot separation is estimated to an accuracy of about V2(0.2), providing 

particle size estimates v^th accuracies on the order of V2(0.2)/A. The error in 
estimating particle size is inversely proportional to the number of pixels separating the 
20 glare spots. Therefore, the CCD camera used for such measurements must have 
sufficient resolution (high pixel count) to distinctly image the glare spots. 

Description of Planar Size and Velocity Measurement Technique 

A technique making use of the characteristics of the glare spots 
predicted by Mie theory is used in the present application of simultaneous 
25 measurement of particle size and velocity using PIV. A PIV system is used to image 
the light scattered from spherical particles at an observation angle of 90**. This 
recorded signal is representative of the Fourier transform of the Mie scattered electric 
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fields incident on the lens of the imaging system. As shown in the previous section, 
the location of the glare spots resulting from reflected and transmitted rays in the 
direction of observation is a good indicator of particle size. Recall that it is crucial to 
have sufficient spatial resolution to image the glare spots. This becomes very 
5 important v^th small particles v^here the glare spot separation is minimal. The 
centroids of the glare spot peaks, and hence the separation between the glare spots, are 
determined using a three-point Gaussian peak estimate technique. Using Equation (6), 
the particle diameter is calculated based on the measured glare spot separation and the 
optical system parameters. Software for performing the size estimations is described 

10 in the next section. Velocity estimates are obtained via PIV correlation processing 
performed on the same images using any commercial or custom PIV processing 
software. Once the size and velocity are obtained, the mass flux of the particle is 
determined according to conventional techniques. Altematively, particle tracking 
velocimetry techniques can be employed, wherein the individual particle velocities are 

15 determined, which can be directly matched to the individual particle size estimates 
obtained fi-om the glare spot separation. 

Sizing 

With reference to FIGURE 5, an algorithm 60 has been developed to 
scan the PIV images for particle images and determine the size of each particle 

20 detected. An image file is read iri a step 62. Thresholding is applied within the image 
(e.g., global thresholding), in a step 64, to reduce noise v^thin the image. Global 
thresholding imposes a user-defined value that limits the minimum non-zero intensity 
value of the image. Any pixels below the threshold value are set to zero (0). The 
thresholding eliminates background noise in the image, thereby making it easier to 

25 discriminate particle images from the background and aiding in the separation of 
overlapping particles. It is important not to pick a threshold level so high that the two 
(2) glare spots are separated from each other. Altematively, a local adaptive 
thresholding technique could also be employed. Local thresholding would offer better 
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image extraction than a global technique, since variations in background levels in the 
image typically require a higher than optimum global threshold level setting. 

The image pixel array is scanned in a step 66. Preferably, the image is 
scanned from left to right, and top to bottom for non-zero (0) pixels. Whenever a non- 
5 zero pixel is found, the eight (8) adjacent pixels are scanned and any non-zero pixels 
are grouped along with the initial pixel to form a particle image array. Then the eight 
(8) adjacent pixels to each of the non-zero pixels are scanned and non-zero pixels are 
grouped into the same particle image array. Pixels cannot belong to more than one (1) 
of the particle image arrays. 

10 Scanning continues until no more adjacent non-zero pixels are found 

and the next particle image is found. Each particle image array is analyzed to 
determine if the image is saturated (i.e., more than one (1) pixel in the array has the 
maximum intensity value as determined by the dynamic range of the camera). The 
sizing algorithm does not accommodate saturated images since it is difficult to 

15 accurately determine the locations of the glare spots v^thin the particle image when 
they are saturated. Therefore, saturated particle images are rejected. 

The glare spot positions and particle diameter are determined in a step 
68. If the particle is not saturated, the maximum intensity in the image is found. The 
values of the pixels above and below, and to the right and left of the maximum 

20 intensity pixel are used to calculate the 3 -point Gaussian estimate of the glare spot 
centroid. Next, the particle image array is scanned to the left and right (or top and 
bottom in the case of a vertical propagating light sheet) to find the point where the 
slope of the intensity distribution curve changes sign, indicating the start of the second 
glare spot. The particle image array is then scanned from the point of slope change to 

25 the edge of the array, continuing in the same scanning direction, in search of the 
maximum intensity within that area of the array. This maximum intensity indicates 
the location of the second glare spot. The same 3 -point Gaussian estimate approach is 
taken to estimate the second glare spot centroid. The difference between the centroid 
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positions represents the variable A in Equation 6 for calculating the particle image 
diameter. This same process is followed for every particle image array within the 
entire PIV image. 

If no slope change is found, then a second glare spot does not exist and 
5 that particle is rejected. The particle may be too small to distinguish between the two 
(2) separate glare spots and, therefore, sizing is not possible using this method. The 
locations of the sized particle images and the calculated diameter are output to a file. 

Particle diagnostics (e.g., diameter, position, distribution mean, and 
standard deviation) are output to a processor in a step 69. 
10 FIGURE 6 illustrates glare spot intensity fields for a range of particle 

diameters. Note that intensity fields 70a, 70b, 72a, 72b, 74a, 74b, 76a, 76b, 80a, 80b, 
82 are very well defined for peirticles having relatively larger diameters. For example, 
the glare spot intensity fields 70a, 70b for a particle having a diameter of 150 microns 
are much more defined than the intensity fields 80a, 80b for a particle having a 
15 diameter of 1 micron. The intensity field 82 only indicates a single glare spot and, 
therefore, would be rejected by the algorithm described above. 

The preferred embodiment of the present invention has been described 
with respect to estimating the particle size as a function of glare spot separations. 
However, it is also contemplated to determine particle sizes as a function of a ratio of 
20 the amplitude of the glare spot intensities. 

Furthermore, it is also contemplated to estimate the diameter of the 
particle by iteratively computing the electric field scattered by a particle and 
determining the resultant electric field detected by the CCD sensor. The process is 
repeated until the correct particle size is determined. Convergence to the correct 
25 particle size is determined when the computed image of the electric field on the CCD 
sensor matches the measured glare spot separation and intensity distribution for a 
given particle image. 
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The invention has been described with reference to the preferred 
embodiment. Obviously, modifications and alterations will occur to others upon 
reading and understanding the preceding detailed description. It is intended that the 
invention be construed as including all such modifications and aherations insofar as 
5 they come within the scope of the appended claims or the equivalents thereof. 
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